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Проанализированы модели расчета высокочастотных вольт-фарадных характеристик. На основе выбранной модели разработа- 
на программа расчета параметров МДП-структуры, а также идеальной вольт-фарадной характеристики и плотности поверх- 
ностных состояний на границе раздела полупроводник-диэлектрик. 


Введение 

Методы исследования структур металл-диэлек- 
трик-полупроводник (МДП) широко развиты, что 
позволяет определить любой параметр, характери- 
зующий качество границы раздела полупроводник- 
диэлектрик и диэлектрического слоя. Наличие на 
границе раздела неоднородностей, обусловленных 
дефектами различного рода, приводит к измене- 
нию энергетического спектра электрона на поверх- 
ности по сравнению с объемом полупроводника. 
Наиболее полная информация о состоянии грани- 
цы раздела содержится в энергетическом спектре 
плотности поверхностных состояний (ППС, 
А ХЕ)). Для его расчета, как правило, используют 
вольт-фарадные характеристики (ВФХ), измерен- 
ные на низкой (НЧ) или ВЫСОКОЙ (ВЧ) частотах. 
Использование для расчета АДЕ) той или иной 
ВФХ имеет свои преимущества и недостатки. 

Основным преимуществом метода НЧ ВФХ по 
сравнению с методом ВЧ ВФХ является возмож- 
ность определения АДЕ) в более широком диапа- 
зоне энергий в запрещенной зоне. В отличие от ВЧ 
метода, метод НЧ позволяет исследовать не только 
в области обеднения-обогащения, но и сильной 
инверсии. Несмотря на это предпочтительнее ис- 
пользовать ВЧ метод, т. к.: 

• Измерения ВЧ ВФХ характеристик технически 
проще. 

• Максимально достижимая чувствительность 
определения АДЕ) в ВЧ методе составляет 
ІО 8 см _2 эВ' ! (в НЧ методе - ІО 10 см' 2 эВ _1 ). 

• Для использования ВЧ свойств полупроводни- 
ков можно использовать характеристики МДП- 
структур с диэлектриком, имеющим сквозные 
токи порядка ІО -9 А/мм 2 (в НЧ методе ток утеч- 
ки не должен превышать ІО 15 А/мм 2 ). 

Существует много методов вычисления плотно- 
сти поверхностных состояний по ВЧ ВФХ (см., на- 
пример [1]). В этой связи необходимо их сравнение 
по выбранным критериям: точности определения 
ППС, анализируемого энергетического диапазона. 

Целью данной работы является разработка про- 
граммы расчета параметров МДП-структуры и по- 
строения зависимости АДЕ) по методу Термана, 
который основан на теоретически измеренной ВЧ 
ВФХ. 


Метод Термана 

В предположении однозарядового поверхност- 
ного состояния (ПС) ППС может быть определена 
по сдвигу экспериментальной ВФХ относительно 
расчетной. Если все зарядовые центры, локализо- 
ванные вблизи границы раздела, как со стороны 
полупроводника, так и со стороны диэлектрика, 
«собрать» на границу раздела, то спектр ППС мо- 
жет быть определен с помощью выражения [1]: 

ЖЛуО = —— Д С) , (1) 

Ч ау* 

где С дюл - емкость диэлектрика, д - элементарный 
кулоновский заряд, у $ - безразмерный поверхност- 
ный потенциал, А Ѵ(С) - разность напряжений 
между теоретической и экспериментальной ВФХ 
при заданном значении емкости. 

Как следует из выражения (1), для построения 
зависимости А, /у,) необходимо выполнить следую- 
щие операции: 

• по экспериментально найденным значениям 
емкости МДП-структуры в режиме обогащения 
и сильной инверсии построить теоретическую 
ВФХ при отсутствии каких-либо ПС; 

• сравнивая экспериментальную и теоретиче- 
скую ВФХ, определить зависимость А К(С); 

• с помощью теоретической ВФХ С(у,) из зависи- 
мости А Ѵ(С) построить А К()Д; 

• продифференцировав по поверхностному по- 
тенциалу зависимость А Ѵ(у 3 ) и используя (1), по- 
строить спектр поверхностных состояний А, /уф 

2. Расчет теоретической высокочастотной ВФХ 

2.1. Классическое решение 

Выражение для теоретической ВЧ ВФХ может 
быть легко получено в рамках модели Гаретта- 
Брэттана исключением из области пространствен- 
ного заряда (ОПЗ) заряда неосновных носителей 
[2]. Это объясняется их неподвижностью в поле ВЧ 
тестового сигнала. Тогда выражение для расчета ВЧ 
ВФХ области пространственного заряда будет 
иметь вид: 

г 5І§п(у д )гг 0 Я' 1 (1-е ѵ -) 

опз 2Ь Д [Г\е у ° -1) + (А-А“ 1 )л] 1/2 ’ 
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где в - диэлектрическая проницаемость полупро- 
водника, - электрическая постоянная, Ь д - дли- 
на Дебая собственного полупроводника, 

А = — = — , где Пі - собственная концентрация 

% п і 

носителей в полупроводнике, р 0 , щ - концентрация 
ионизованных акцепторов, доноров в примесном 
полупроводнике. 

Рассчитанный по (2) график зависимости ВЧ 
ВФХ от показан на рис. 1 , а. Выражение (2) полу- 
чено в предположении, что подвижные носители 
заряда подчиняются статистике Максвелла-Боль- 
цмана, но т. к. при больших поверхностных потен- 
циалах^ возможно их вырождение, то это приведет 
к ошибке в определении С опз , рассчитанной по (2). 
Следовательно, в указанных областях поверхност- 
ного потенциала (обогащение, инверсия - сильная 
инверсия) величина расчетной емкости будет завы- 
шена по отношению к реальной. Это приведет к 
ошибке в определении ППС. 


2.2. ВФХ характеристика с использованием 
распределения Ферми-Дирака 

Согласно [3] выражение ВЧ ВФХ для полупро- 
водника п типа на участке обогащения и обеднения 
выглядит следующим образом: 


С от (Л ) = 4§п( ѵ ) 
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где 8 й — отношение краткости вырождения основ- 
ного состояния ионного остатка к краткосрочно- 
сти вырождения основного состояния нейтрально- 
го донора, Р у / 2 ( 77 ), Ру 2 (гі) ~ интеграл Ферми- Дирака 
степени 1/2 и 3/2, \Ѵ 0Р — энергия от уровня энергии 
донора до уровня Ферми в легированном полупро- 
воднике, ]Ѵ К и \Ѵ, С - энергия от уровня Ферми в ле- 
гированном полупроводнике и энергия от уровня 
Ферми в собственном полупроводнике до дна зоны 
проводимости, у іпѵ - соответствует границе между 
инверсий и обеднением. 

Расчет ВЧ ВФХ по (3) представлен на рис. 1, б. 
Видно, что при использовании распределения 
Ферми-Дирака при больших^ ВФХ становится бо- 
лее плавной по сравнению с (2). Изменение стати- 
стики при больших у \ по очевидным причинам 
приводит также к большему значению емкости 
ОПЗ в режиме «инверсия - сильная инверсия». 
Недостатком расчета ВЧ ВФХ этим методом явля- 
ется необходимость вычисления огромного коли- 
чества интегралов Ферми-Дирака, например, с по- 
мощью полиномов Лагерра. 


Согласно [2], максимальная ошибка в опреде- 
лении емкости ОПЗ при 300 К для минимального 
значения емкости ОПЗ равна 100 %, а для емкости 
плоских зон - 23,5 %. 


2.3. ВФХ кривая с перераспределением носителей 
заряда инверсионного слоя 


Данная модель рассчитывается при следующих 
физических предположениях [4]: 

• Равновесие по постоянному смещению. Для из- 
мерений ВФХ это эквивалентно тому, что ско- 
рость нарастания напряжения недостаточна для 
избежания глубокого обеднения. 

• Учет внутреннего неравновесия носителей за- 
ряда в инверсионном слое, связанного с их вы- 
сокой инерционностью и малым периодом те- 
стового ВЧ сигнала. Это приводит к эффекту 
перераспределения носителей заряда, который 
не учитывался в предыдущих расчетах. 

Согласно [4] выражение для ВЧ емкости обла- 
сти пространственного заряда в полупроводнике р 
типа имеет вид: 

г 

^ ОПЗ 
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где С ш 7 - низкочастотная емкость ОПЗ в полупро- 
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Расчет ВФХ по (4) представлен на рис. 1, в . Полу- 
ченное значение емкости хорошо совпадает во всем 
диапазоне с методом, представленным в п. 2 . 1 . 
Принципиальным недостатком этого метода расчета 
является наличие глубокого минимума вблизи точки 
(функция меняет свой знак вблизи указанной точки) 
(рис. 1,в). Использование данного метода для расче- 
та ППС без операций сглаживания кривой невоз- 
можно ввиду принципиального искажения ВФХ. 

Таким образом, для вычисления спектра ПС в 
режиме реального времени наилучшим образом 
подходит расчет ВЧ ВФХ по выражению (2). Воз- 
никающая в результате его использования ошибка 
в спектре ПС существенна лишь в области обога- 
щения и сильной инверсии, где ППС, по методу 
Термана, рассчитывается с большой ошибкой [5]. 


3. Описание программы 

Используя методику расчета ППС МДП-структу- 
ры, написана программа, позволяющая на основе экс- 
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Рис. 1 . ВЧ ВФХ зависимости ОПЗ в идеальной МДП -структуре: а) классическая, б) с применением распределения Ферми-Ди- 
рака, в) с учетом перераспределения носителей заряда инверсного слоя 


периментальных ВЧ ВФХ определять N ХЕ), а также 
основные параметры МДП-структуры. Вид пользова- 
тельского интерфейса программы показан на рис. 2. 

Окно программы имеет четыре вкладки: 

Вкладка «Измерение» 

На ней находится набор элементов для управле- 
ния измерительным комплексом «МЕТРО- 
НОМ-2», описанным в работе [6]; 

Вкладка «Расчет параметров МДП-структуры» 

Здесь осуществляется открытие ранее сохранен- 
ных данных при нажатии кнопки «Открыть». Име- 
ется возможность одновременного построения 10- 
ти графиков. Если открытые данные ВФХ не 
устраивают оператора по каким-либо причинам, то 
их можно обнулить, указав номер обнуляемого гра- 
фика в элементе СотЬоВох с названием «График» и 
нажав кнопку «Сброс». Построенные графики мож- 
но сохранить в виде массива чисел. Так же возмож- 
но распечатать отображаемые результаты ВФХ. 

Перед началом расчета параметров МДП-струк- 
туры необходимо указать следующее: материал по- 
лупроводника, материал диэлектрика, площадь по- 
левого электрода МДП-структуры, температура, 
при которой проводилось измерение. Путь к по- 
следнему открытому файлу с данными отображает- 
ся в поле «Путь». Расчет параметров МДП-структу- 
ры осуществляется нажатием кнопки «Расчет». 

Рассчитывается выбранный график согласно 
алгоритму, изложенному выше. Для того чтобы из 


расчета исключить случайные ошибки в величине 
емкости перед началом расчета параметров МДП- 
структуры осуществляется согласование экспери- 
ментальной ВФХ при помощи алгоритмов сглажи- 
вания. Чувствительность сглаживающего алгорит- 
ма к различного рода ошибкам на ВФХ задается в 
текстовом поле с названием «Коэффициент сгла- 
живания». На рис. 2 показана экспериментальная 
кривая, полученная при коэффициенте сглажива- 
ния ІО -3 и отмеченная численными значениями 
С диэл , емкости плоских зон и емкости в инверсии. 

Из сглаженной экспериментальной кривой ВЧ 
ВФХ определяется емкость диэлектрика, емкость 
плоских зон, а также минимальная емкость ОПЗ и со- 
ответствующие им напряжения. Полученные значе- 
ния емкостей отображаются на сглаженной ВФХ кри- 
вой в виде значений емкостей (рис. 2). Кроме этого, 
используя данные о МДП-структуре, рассчитываются 
следующие параметров: тип проводимости полупро- 
водника, толщина диэлектрика, максимальная шири- 
на ОПЗ, концентрация доноров (акцепторов), ем- 
кость плоских зон, напряжение плоских зон, встроен- 
ный заряд в диэлектрике и на поверхностных состоя- 
ниях границы раздела полупроводник-диэлектрик. 

Используя рассчитанные параметры, строится 
теоретическая ВФХ по формуле (2). Как показано 
выше, данное выражение оптимальным образом 
описывает емкость ОПЗ в идеальной МДП-струк- 
туре. Пересчитанная функция емкости от поверх- 
ностного потенциала С(у 5 ) в функцию от напряже- 
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Рис. 2. Вид пользовательского интерфейса 

теоретическая и экспериментальная С-Ѵ у(Ѵ) 




сіеПаѴ(Ѵ) 



Рис. 3. Вкладка «Промежуточные результаты» 


ния С(Ѵ) отображается на графике вместе с экспе- 
риментальными данными. 

Вкладка «Промежуточные результаты» 

Для контроля правильности расчета параметров 
на вкладке «Промежуточные результаты» выводят- 
ся три графика (рис. 3): теоретическая и экспери- 
ментальная сглаженная ВФХ, зависимость поверх- 


ностного потенциала у 5 от приложенного напряже- 
ния, разность напряжений между теоретической и 
экспериментальной кривыми при заданном значе- 
нии емкости от приложенного напряжения. 

Вкладка «График 7Ѵ„» 

Здесь (рис. 4) отображается рассчитанный по 
ур. (1) спектр ППС МДП-структуры. 
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Зависимость N53 [1 /эВ см л 2] от Е [эВ] 



Рис. 4. Вкладка «Г рафик Л/ 55 » 


Для качественной оценки рассчитанных дан- 
ных на графике отмечены энергии, соответствую- 
щие дну зоны проводимости (Е с ), потолку валент- 
ной зоны ( Е ѵ ), уровень Ферми (Е), а также симме- 
тричному уровню Ферми относительно такового в 
собственном полупроводнике (2 Е г ), уровню Ферми 
в собственном полупроводнике (Д), уровню энер- 
гии, при которой наступает глубокая инверсия 
(21п(А)-3). Полученная зависимость 7ѴДД) может 
быть сохранена в виде точечного рисунка и распеча- 
тана на принтере. Как можно видеть, полученный 
спектр ППС имеет классическую [/-образную фор- 
му и типичные значения ППС для кремниевых 
МДП-структур. В дополнение к этому рассчитан- 
ный спектр ППС демонстрирует наличие локальных 
минимумов ППС вблизи середины запрещенной 
зоны полупроводника. Данное обстоятельство явля- 
ется спецификой использованной МДП-структуры, 
связанной с особенностями технологии ее изгото- 
вления и обработки. 


Заключение 

Разработана программа расчета параметров 
МДП-структуры, позволяющая на основе экспери- 
ментальных высокочастотных ВФХ определить 
тип проводимости полупроводника, толщину ди- 
электрика, максимальную ширину ОПЗ, концен- 
трацию доноров (акцепторов), емкость и напряже- 
ние плоских зон, величину суммы встроенного за- 
ряда в диэлектрике и на поверхностных состояниях 
раздела полупроводник- диэлектрик. 

Используя параметры МДП-структуры, про- 
грамма позволяет рассчитать теоретическую ВФХ 
и спектр плотности поверхностных состояний, что 
может иметь самостоятельное применение. 

Все результаты расчетов отображаются в графи- 
ческой форме и могут быть сохранены или распе- 
чатаны на принтере. Программа имеет удобный 
пользовательский интерфейс с возможностью 
просмотра промежуточных результатов расчета. 
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